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Based on DC magnetization measurements, the temperature dependencies of the trapped magnetic 
field have been calculated for two MgB2 samples prepared by two different techniques: the high- 
pressure sintering and the hot pressing. Experimentally measured trapped field values for the first 
sample coincide remarkably well with calculated ones in the whole temperature range. This proves, 
from one side, the validity of the introduced calculation approach, and demonstrates, from another 
side, the great prospects of the hot pressing technology for large scale superconducting applications 
of the MgB 2 . 

PACS numbers: 74.25.Sv, 74.70.Ad 
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Magnesium diboride is a new promising superconduct- 
ing material with a critical temperature of about 40 K. 
High T c values and simple chemical composition of MgB2 
have made it an interesting object for applied investiga- 
tions with great perspectives of use in superconducting 
motors, flywheels and bearings. The key parameter for 
such applications is the maximum trapped field in the 
sample and its temperature dependence. This, in turn, 
is closely connected both with the critical current density 
and with the size of superconductor. 

Since the very discovery of superconductivity in MgB2 
the researchers dealt with tiny samples which were not 
suitable for large scale applications. Now, the situation is 
changed and several techniques have been elaborated to 
produce high quality MgB2 bulk polycrystalline samples 
of a few centimeters in diameter and with the critical 
current densities up to 10 5 - 10 6 A/cm 2 0,0- 

It is remarkable that, in contrast to the high tempera- 
ture superconductors, the grain boundaries in bulk mag- 
nesium diboride superconductors do not act as weak links 
3, 4J. This significantly simplifies the growth of bulk 
samples suitable for large scale applications. 

In this letter we focus on an important property of the 
bulk superconductor - the trapped magnetic field. We 
calculate its temperature dependencies from DC magne- 
tization data using an iteration approach. We show that 
the trapped field values measured experimentally are in a 
perfect agreement with the calculated ones. This makes 
us able to predict the expected values of the trapped 
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FIG. 1: Temperature dependencies of maximum trapped field 
in bulk MgB 2 samples prepared under high-pressure sintering 
and ball-milled techniques. 



magnetic field in a sample which is not yet grown to an 
appropriate size but which preparation technique looks 
the most promising today. 

We have studied two MgB2 bulk samples prepared by 
different techniques. One of them was sintered under 
high pressure as described in p|. The sample has an 
uniform polycrystalline structure without cracks and has 
28 mm in diameter and 11 mm in height. The trapped 
magnetic field in the centre of the sintered sample was 
measured by a Hall probe for different temperatures from 
6 to 33 K. The resulting temperature dependence of the 
trapped field has a negative curvature and trends to sat- 
uration at low temperatures (see Fig. 1, cross symbols). 
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FIG. 2: Field dependencies of critical current density at dif- 
ferent temperatures for high- pressure sintered MgB2 sample. 
J C {H) curves were obtained from DC magnetization data. 

The other sample was prepared by ball milling of Mg 
and B powders at ambient temperatures followed by hot 
pressing [(|, which we call " ball- milled" . It consists of 
spherical nanocrystalline grains about 40 - 100 nm in size, 
that distinctly improve pinning due to the large number 
of grain boundaries. For this sample the trapped field 
can not be measured because of its small size. 

To calculate the trapped field one should know tem- 
perature and field dependencies of the critical current 
density J C (H,T) of the both superconductors. For DC 
magnetization measurements small bar shaped pieces of 
the both samples have been used. Then, the J C {H) curves 
were obtained from the magnetization loops M(H) using 
the conventional expression j c (H) — 20 • AM(b — b 2 /3l), 
where AM is the difference of the magnetization (in 
emu/cm 2 ) measured for ascending and descending ap- 
plied field, b and / are the sample width and length in 
cm, respectively, and j c is obtained in A/cm 2 . 

For the high-pressure sintered sample a set of J C {H) 
curves presented in Fig. 2 in logarithmical scale 
are straight lines for low fields up to the level of 
J c = 10 4 A/cm 2 . This means an exponential decay 
of critical current density with increasing field 

J C (H, T) = J(T)exp{-H/b(T)), (1) 

where J(T) and b(T) are temperature dependent coeffi- 
cients. 

The dependencies of J(T) and b(T) for the sintered 
sample, shown in Fig. 3 by circles, can be well fitted to 
the expressions 

b(T) = 6 [1 - {T/T c )% J(T) = J [l - (T/T c )]. (2) 

Here T c = 37.5 K is the critical temperature, whereas 
b = 0.9 T and J = 3.1 x 10 5 A/cm 2 are fitting parame- 
ters, which characterize an effective field of supercurrent 
decay and a maximum critical current density at zero 
field and temperature, respectively. 
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FIG. 3: J(T) and b(T) coefficients of J C (H) curves (see Eq.QJ 
for MgB2 bulks. Symbols represent experimental data for 
high-pressure sintered (circles) and ball-milled (squares) sam- 
ples. Solid lines represent fitting results, which is described 
by formulas (5J and for sintered and ball-milled samples 
respectively. 



For the ball-milled MgB2 sample J C (H,T) can be de- 
scribed by Eq. Q too (see Fig. 4 in but the functions 
b(T) and J(T) are different: 

b(T) =bo[l- (T/T c ) 2 ] 3 / 2 , J(T) = J [l - (T/T c ) 2 ] 3 / 2 , 

where the corresponding parameters are T c — 34 K, 
b Q = 2.78 T and J = 8.5 x 10 5 A/cm 2 (see Fig. 3, 
squares). 

The obtained results show that the critical current 
density strongly depends on temperature. For the high- 
pressure sintered sample the values of J(T) and b(T) are 
about 2.5 times lower than those of for the ball- milled 
sample. From the other hand, the latter sample has a 
much stronger power dependence of J(T) and b(T) on 
temperature. 

The further calculation of the trapped field from a 
known J C (H,T) function is a nontrivial task. According 
to the Biot-Savart low, the magnetic field H generated 
by the supercurrent J c flowing in the volume of sample 
V is 

H(r) = - f J c (H,p) x , r ~ P ,^ p. (4) 
cJv \ r - P\ 

As far function of H, one needs to solve this 

integral equation in order to calculate the trapped field 
values. 

To simplify the problem we assume the homogeneity 
of the sample, i.e. (1) that the critical current density 
does not depend on coordinate p explicitly and (2) that 
the currents flow along concentric circles. Then, the 
field component normal to the sample surface, H z , can 
be written in cylindrical coordinates r — (r, <f>, z) and 
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(p, ip, C) as follows 



where f(r,p,tp,C) 



Jc(H z )f(r,p^,()pd 3 p, 



(5) 



p— r cos ip 



[(r-pcosi/>) 2 + (psini/>) 2 + (z-C) 2 ] 3/2 ' 

d 3 p = dpdipd^. The resulting field profile has a cylin- 
drically symmetric form with the maximum at r = 0. 

The integral equation (JSJ can be solved for i/(r) nu- 
merically using an iterative procedure 



H i+1 (r) = - [ J c (H l )f(r,p^,Opd 3 p, 
c J v 



(6) 



until — Hi\/Hi < e . Here the obtained profile is 

used to calculate current in the next step of iterations. 
We chose an initial profile Ho(r) = const to start the 
iterations. 

One should note that we calculate the trapped field 
profile numerically in each step and can not define H(r) 
function in every point. Therefore, we should discretize 
Eq. JSJ) and divide the cylindrical sample into N concen- 
tric tubes of width A = jrti and then calculate the field 
in each point rk, k = . . . N — 1 as a sum of fields from 
all tubes 

A N-1 

H ( r k) = -J2 MB0>j))f(rk,Pj,1>,Qpjdil><K, (?) 
° 3=0 J 

where rk = kA, R s is the radius, L is the height of the 
cylindrical sample and the integration here is made over 
ijj from to 2-7T and over £ from —L to 0. 

From equations J7J one obtains the field profile for a 
given J C (H) function using iteration (JSJ. In the case of 



J C (H) = const we will obtain a conical Bean profile with 
the maximum in the center of the sample. To reach the 
accuracy of e = 0.1 % it is sufficient to make just 16 
iterations with N = 14 tubes. Moreover, the method 
of calculation has no fitting parameter and use only ex- 
perimental J C {H,T) data. The calculated temperature 
dependence of maximum trapped field for high-pressure 
sintered and ball-milled MgB2 samples are represented 
on Fig. 1 by the solid and dashed lines respectively. 

We obtained a very good coincidence with experimen- 
tal points for the high-pressure sintered sample. Remark- 
ably, the trapped field, which is determined by the cur- 
rent distribution in the whole sample volume can be cor- 
rectly calculated on the base of J c measured locally in 
a small piece of the sample. The ball-milled (nanocrys- 
talline) sample was found to have much higher trapped 
fields than the sintered sample. Thus, the development 
of large size nanocrystalline MgB2 samples is very attrac- 
tive for future applications. 



To sum up, we calculated the temperature dependence 
of the trapped magnetic field for two MgB2 bulk samples 
prepared under different techniques. We used DC mag- 
netization data and an iteration approach to solve the 
Biot-Savart equation. Experimentally measured trapped 
field values for the high-pressure sintered sample coincide 
excellent with calculated ones in the whole temperature 
range. The correctness of the calculation allowed us to 
predict the expected values of the trapped field in the 
ball-milled sample and revealed the great prospects of 
this technology for large scale applications of supercon- 
ducting MgB2. 
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